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A  THEORETICAL  INVESTIGATION  OF 
HYDROGEN- AIR  REACTIONS 

PART  I  -  BEHAVIOR  WITH  ELABORATE  CHEMISTRY* 

By  P0  A.  Libby1,  H.  Pergament*  andM.  H.  Bloom1 

SUMMARY  .  . 

As  part  of  a  study  of  hydrogen-air  reactions  under  conditions  of  interest  in 
propulsive  devices,  a  theoretical  analysis  of  these  reactions  with  an  elaborate 
description  of  the  chemical  behavior  considered  has  been  carried  out,  The 
system  has  been  idealized  so  as  to  result  in  an  adiabatic,  isobarlc  process. 

In  the  future,  models  with  simplified  chemistry  will  be  treated  in  order  to 
establish,  if  possible,  a  chemical  kinetic  description  which  can  be  employed, 
with  greater  ease  than  the  chemical  reactions  considered  in  this  report,  for  the 
treatment  of  practically  interesting  flows.  The  present  analysis  yields  the 
temperature  and  composition  histories  for  a  range  of  initial  temperatures, 
mixture  ratios  and  pressures. 


*  This  research  was  sponsored  by  the  U»  S.  Air  Force  through  the  Office  of 
Scientific  Research,  Contract  No.  AF  49{638)-99l.  Technical  monitoring 
was  done  by  Dr.  Milton  M.  Slawsky,  Chief,  Propulsion  Research  Division, 
AFOSR* 

1  Senior  Scientific  Investigator 

2  Senior  Scientist 
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SYMBOLS 

constants  defined  by  Equation  {14} 

1 

B 

! 

0 

concentration  of  mixturet  mol/ee 

j 

ci 

concentration  of  species  i*  mol/cc 

j 

CP 

specific  heat  at  constant  pressure 

i 

h 

enthalpy  of  mixture  .(defined  by  Equation  (i>)) 

,  ■ 

hi 

enthalpy  of  species  1 

i 

k 

reaction  rate  parameter 

S 

equilibrium  constant 

i 

i 

\ 

P 

refers  to  chemical  symbol  for  products  in  a  reaction 

\ 

P 

static  pressure 

’  i 

i 
i 

R 

refers  to  chemical  symbol  for  reactants  in  a  reaction 

i 

i 

i 

Ru 

universal  gas  constant 

f 

1 

t 

I 

i 

t 

time 

i 

t 

T 

static  temperature 

) 

i 

W 

molecular  weight  of  mixture  {defined  by  Equation  {?)} 

i 

W| 

molecular  weight  of  species  i 

! 

> 

*1 

mass  rate  of  formation  of  species  i 

|- 

X 

mole  fraction 

T 

mass  fraction 

! 

n 

molar  rate  of  formation  of  species,  i 

r 

enthalpy  at  reference  temperature  of  species  i 

r 
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n  equivalence  (mixture)  ratio 
V*  stoichiometric  coefficient  of  react  ant 

Vn  stoichiometric  coefficient  of  product 

p  density  of  mixture 

Tj  denotes  end  of  induction  period  of  reaction 

T  denotes  end  of  second  period  of  reaction  (see  discussion  of  results) 

denotes  time  at  which  equilibrium  is  achieved 
rd  denotes  ignition  delay  time  (see  discussion  of  results  and  Ref.  10) 
\l  sec.  microseconds 
Subscripts 

i  pertaining  to  species  i 

j  pertaining  to  species  j 

m  pertaining  to  reaction  m 

0  condition  at  time  t  c  0 

r  reference  value 


r  :■ 


JL:-. 


J  . , 


v> 
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The  chemical  kinetics  of  hydrogen-air  and  hydrogen-oxygen  reactions  have 
been  studied  extensively  in  the  past.  Reference  1  provides  extensive  information 
on  these  reactions  while  References  2  and  3  are  considered  to  provide  the  most 
recent  and  accurate  data  on  the  mechanism  and  rate  constants  for  gas  phase 
reactions  between  hydrogen  and  oxygen  at  relatively  high  temperatures,  Most  of 
the  studies  of  this  chemical  system  have  been  performed  in  connection  with  constant 
volume  vessels  and  have  involved  significant  heterogeneous  reactions. 

There  has  recently  been  considerable  interest  in  this  reaction  on  the  part  of 
aeronautical  scientists  because  of  the  possible  use  of  hydrogen  as  a  fuel  in  high 
performance  propulsive  devices  (c,£,  References  4-6).  These  applications  involve 
flow  systems  in  which  the  chemical  and  fluid  mechanical  aspects  are  intimately 
related.  The  limiting  models  of  chemical  behavior,  i.e,,  frozen  and  equilibrium 
flows,  are  usually  studied  first  in  order  to  minimize  the  coupling  between  the  two 
aspects.  Ultimately,  however,  the  selection  of  the  more  realistic  of  these  solutions 
requires  simultaneous  consideration  of  the  chemical  kinetics  and  fluid  mechanics. 
The  process  of  the  formation  of  water  from  hydrogen  and  oxygen  has  been 
found  to  be  complex  and  to  involve  many  elementary  reactions.  Indeed,  Duff  and 
Schott  {References  2  and  3)  list  the  following  gas  phase  reactions? 

+  O.  2QH 

H+Oj  yst  OH  4-  0 

O  4  H2  OH  +  If 

OH  +H,==iK10  H*H 
% s  a 

20H  ZS-^H.O  4  O 


•r' 


i  ,  • 


m 


f 
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H  +  M 
2 


H  O  +  M 
2 


(2) 


2H  +M: 

H  +  OH  +  Mj 
H  +  O  +  +  M 

20  +  +  M 

Equations  (l)  proceeding  to  the  right  describe  the  bimoiecular  reactions  dominating 
the  early  phase  of  the  reaction  and  Equations  (2)  proceeding  to  the  right  describe  the 
recombination  reactions  which  dominate  the  final  approach  to  equilibrium,  The 
symbol  M  in  Equations  (2)  denotes  a  third  body* 

In  the  usual  fluid  mechanical  treatment  of  the  limiting  cases  of  the  hydrogen-air 
reaction  in  the  temperature  range  from  1000  to  3000°K,  it  is  customary  to  assume 
as  species  present  in  significant  amounts  02,  Hg  ,  H2  O  and  N2  .  Consideration  of 
Equations  <1)  and  (2)  indicates  that  even  if  nitrogen  is  considered  an  inert  dilutent, 
there  are  three  additional  species,  O,  H,  and  OH  which  must  be  considered  in  the 
study  of  non- equilibrium  flows. 

In  view  of  this  complexity  it  is  important,  for  aeronautical  applications,  to 

consider  simplified  but  reasonably  accurate  descriptions  of  the  chemical  kinetics 

of  the  hydrogen-oxygen  reaction.  From  observations  of  hydrogen-air  flows  Fine 

(Reference  7)  deduced  an  overall  reaction  rate  equation  which  can  be  interpreted  as 

(3/2)  H2  +C2e=3H2  O  +OH 
H2  -f-02p^2OH 

These  measurements  were  roughly  ia  the  range  of  temperatures  between  1300  and 
200 0°K  and  thus  correspond  to  temperatures  of  interest  in  some  aeronautical 


applications  (c£«  Reference  6),  The  appearance  of  the  (3/2)  factor  in  Equation  (3) 
indicates,  according  to  the  usual  chemical  usage,  that  this  equation  provides  an 
overall  description  of  the  production  of  water. 


If  an  overall  reaction  equation  of  the  type  given  hj  Equation  (3)  is  found 
sufficiently  accurate,  then  only  one  additional  species,  OK,  and  two  r ©action 
equations  need  be  considered  in  the  treatment  ©f  Ida®  non- equilibrium  fluid 
mechanics,  and  coupling  between  flow  and  chemical  effects  can  then  be  taken 
into  account  with  a  reasonable  amount  of  effort. 

In  order  to  provide  results  which  can  be  used  to  evaluate  the  accuracy  of 
either  Equation  (3)  or  similar  approximate  descriptions  of  the  reaction,  a 
theoretical  investigation  of  a  model  problem  with  the  elaborate  chemistry  of 
Equations  (1)  and  (2)  has  been  carried  out  and  is  described  here.  In  subsequent 
reports  the  comparison  between  the  results  presented  here  and  the  predictions  of 
approximate  representations  will  be  given.  The  fluid  mechanical  model  considered 
involves  a  constant  pressure,  adiabatic  system,  A  homogeneous  mixture  of  hydrogen, 
and  air  at  specified  values  of  the  initial  temperature,  pressure  and  equivalence 
ratio,  which  do  not  correspond  to  an  equilibrium  ©tat®,  is  assumed  to  prevail  at 


time  zero.  The  subsequent  history  of  the  composition  and  temperature  at  constant 
pressure  is  studied.  Note  that  time  "zero"  can  be  associated  with  a  sudden 
chemically  frozen  compression  ox  expansion  of  an  equilibrium  mixture  to  the 
conditions  corresponding  to  the  initial  state,  j 

i' 

This  model  has  several  physical  interpretation#:  ©n@  correspond®  to  a 
chemically  inert  balloon  surrounding  a  homogeneous  gas  mixture,’  the  second,  to 
a  constant  pressure  expansion  in  a  nozzle  in  which  the  flow  is  treated  a#  one- 
dimensional  and  inviscid;  and  a  third,  to  the  impulsive  motion  of  a  piston  fa  a 
containing  cylinder  in  which  the  pressure  and  initial  temperature  adjustments  j 

occur  prior  to  the  initiation  of  the  chemical  changes.  If  the  pressure  is  I 


maintained  constant  during  the  reaction  in  these  cases,  the  mass-velocity  is  zero 
ore  uniform  and  may  be  disregarded  in  the  energy  balance.  If  desired,  the  area 
distribution  of  the  nozzle  or  streamtubes  can  be  computed  a  posteriori.  It  is  of 
imterest  to  note  that  these  physical  interpretations  permit  the  volumetric  require¬ 
ments  for  the  combustion  of  premixed  hydrogen-air  reactants  to  be  estimated 
paro'vided  &  constant  velocity  prevails. 

This  model  was  selected  for  theoretical  investigation  since  it  corresponds 
miore  closely  to  the  conditions  of  interest  in  connection  with  diffusive  burning  of 
hydrogen  (cf.  Reference  6)  than  the  case  usually  treated  by  chemical  kineticists. 
Moreover,  the  model  is  similar  to  problems  studied  by  Duff  (Reference  2)  who 
determined  the  reaction  profile  behind  detonation  waves  in  hydrogen-oxygen  mixtures 
diluted  in  xenon.  However,  such  reaction  zones  involve  significant  changes  in  both 
velocity  and  pressure  during  the  period  of  reaction.  It  should  be  noted  that  if  the 
initial  state  is  assumed  to  be  achieved  by  a  strong  shock,  the  pressure  alterations 
do-vnstraam  of  the  shock  will  be  small  and  the  model  employed  here  will  be 
ap»pro3cimately  valid.  This  shock  in  general  will  not  satisfy  the  Chapman-Jouguet  con 
di trio n  but  can  be  created,  for  example,  by  the  high  speed  motion  of  a  blunt  body 
thorough  &  hydrogen-air  mixture. 

In  this  analysis  it  is  assumed  that  the  maximum  temperatures  are  sufficiently 
low  so  that  nitrogen  may  be  considered  an  inert  dilutent,  i»e=,  the  dissociation  of 
nitrogen  and  the  formation  of  nitric  oxide  and  ammonia  are  neglected.  Thus  the 
sp aerie s  considered  are  02»  H3,  H20,  N2,  0,  H  and  OH;  these  will  be  denoted 
respectively  by  subscripts  Is 

In  the  following  section  the  theoretical  analysis  is  described  while  in  the  final 
sections  the  numerical  results  and  conclusions  are  discussed. 
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Auxiliary  equations  of  the  form 


7  7 


\  VI,  \  *w4-o 

m 

^  =  Vi^  Oj ,  T),  iijsli  2,  ...  7 

(9) 

hj  =>  tyT) 

uo> 

Cj-PYj/Wj 

(id 

complete  the  system  of  equations.  The  initial  conditions  are:  at  t  *  0, 


.  'j-  's'  **  '  ■ ' 

.  . .  _  !  f 


■  v 


Consider  the  detailed  expressions  for  the  y.’s;  first,  with  nitrogen 
assumed  to  be  inert,  y4  =  0  and  Y4  =  Next,  since  the  activation  energy 


of  the  first  reaction  equation  in  the  system  given  by  Equation  (1)  is  consider* 


ably  larger  than  the  others,  it  is  reasonable  to  neglect  this  reaction  and  to 


consider  only  the  last  eight  reactions.  The  rate  constants  for  the  forward 


reaction,  i.e. ,  in  the  direction  to  the  right  in  the  Equations  (1)  and  (2),  are 


denoted  km,  m  =  1,  2. .  .8  while  the  corresponding  backward  constants  are 


denoted  km,  m  =  9,  10, ...  16.  The  general  expression  for  Vj  is  then  (c£« 


Reference  8) 


(V  "  -  V\  )  k 

i»m  i,  m  rn 


.  V\ 

•  -XA 

Ci  ,  l,  j#4 


for  the  general  reaction  ^  ^  V'.' 


The  values  of  the  rate  constants  for  the  forward  reactions  were  taken 


from  Schott  (Reference  3);  those  for  the  backward  reactions  were  computed 


from  the  equilibrium  constants  corresponding  to  each  elementary  reaction* 


Thus  (cf.  Reference  8) 


i  An, 

k  =  k  K”1  (R,  T)  m 

m  m-8  P,m  u 


where  m  =  9,  10, ...  16  and  An, m= 


(V"  -  V\  ) 

i,  m  i,  m 


The  equilibrium  constant  K  in  Equation  (13)  was  approximated  la 

P*  33*1 


form 


“X- 


i  V 


!  *  •  V 


E 


Am/T 

^p,  m  s  8 

over  the  temperature  range  1GQQ°K  <  T  <  3000°K;  the  requisite  values  of  K^m 
being  computed  by  application  of  the  tables  in  Reference  8.  In  Table  1  and 
Appendix  I  there  are  presented,  respectively,  the  values  for  k^  and  the  explicit  ‘ 

relations  for  the  y^’s  used  for  this  analysis.  The  third  bodies  in  the  recombina-  J 

tion  reactions  given  by  Equations  (2)  were  taken  to  be  all  species  present.  f 

Consider  next  the  determination  of  the  initial  conditions.  It  is  assumed  that  ’ 

j 

| 

| 

at  time  zero  the  atomic  hydrogen  and  atomic  oxygen  are  in  equilibrium  with  H2  and  02  j 

•1 

respectively.  Moreover,  neither  water  nor  hydroxyl  radical  is  assumed  present 
so  that  at  t  =  0,  Y  ^  =  Yj  =  0.  The  initial  concentrations  of  molecular  oxygen  and 
molecular  hydrogen  and  the  constant  concentration  of  nitrogen  are  conveniently 
expressed  in  terms  of  the  equivalence  ratio  7}  defined  by  the  chemical  equation 

02  +  3.76N2+  lr)  H2  products  (15) 

Thus  a  value  of  r;  -  1  corresponds  to  stoichiometric  proportions  of  hydrogen  and 
oxygen.  In  terms  of  7]  the  remaining  initial  conditions  are  j 

Y1(  0  =  32/(137.  3  +  4t}) 

y2,0  =  477/(137.3  +  41?)  (16) 

Y4}0  =  105.  3/(137,3  +  4f?) 

Clearly  this  does  not  imply  equilibrium  of  the  entire  mixture  in  the  initial 
state. 


Now  Equation  {16}  along  with  the  stated  initial  conditions  on  and  Y^ 
permit  the  approximate  initial  concentrations,  Y^q  and  g  to  be  obtained 
from  the  equilibrium  constants  for  the  dissociation  of  molecular  oxygen  and 
hydrogen  respectively. 

It  is  noted  that  other  treatments  of  the  initial  conditions  could  have  been 
employed.  For  example,  it  could  have  been  assumed  that  Y^  =  Y,.  =  Y^  =  Y  =  0 
at  t  =  0.  Since  the  first  reaction  in  Equation  (1)  was  neglected,  the  reaction 
would,  in  this  case,  be  initiated  by  dissociation  of  H2  and  02  followed  by  the 
chain  reactions  leading  to  the  formation  of  OH,  and  additional  O  and  H.  It  is 
believed  that  this  treatment  would  alter  only  the  induction  period  described 
below  and  would  not  affect  significantly  the  overall  reaction  history. 

Numerical  analysis  requires  analytic  representation  of  the  static 
enthalpy-temperature  relations  corresponding  to  Equation  (10).  For  the 
limited  temperature  range  of  interest  in  this  study,  it  will  be  sufficiently 
accurate  to  let 


^  =  &i  +  5p.(T  -  Tr) 

where  Tr  is  a  reference  temperature  in  the  range  of  interest  and  and 
are  constants  so  that  the  exact  =  h^(t)  relation  is  closely  approximated. 
The  assumption  attendant  with  Equation  (10a)  implies  rigorously  a  frozen 
vibrational  energy  at  some  mean  temperature.  Equation  (10a)  can  be 
employed  in  Equation  (6)  to  yield  the  following  explicit  expression  for  the 
temperature: 


i 


(10a) 


12 


■  n 


,  *  i 


y.  V 


T  ■  Ty  +  j 


Yi  h^pJT-Ty)  j; 


Th®  analvsi a  required  for  numerical  determination  of  th©  history  of  a 


mixture  of  hydrogen  and  air  according  to  the  model  under  consideration  has 


now  been  presented.  A  particular  problem  is  specified  by  value®  of  the 


pressure  p,  th®  initial  temperature  TQ>  and  th®  equivalence  ratio  77 


The  numerical  integration  was  carried  out  on  an  IBM  7090  employing 


the  standard  fourth  order  Runge-Kutta  method  for  starting  the  calculations 


and  the  Adam  s  -  Mo  ul  to  n  predictor-corrector  technique  for  the  continuinf 


integration.  The  solution  of  the  differential  equations  was  obtained  to  times 
of  roughly  2(10“^)  seconds  or  more,  or  until  the  temperature  became  essentially 


constant  with  time. 


It  is  noted  that  while  in  principle  the  numerical  integration  may  be 


carried  out  in  terms  of  five  differential  equations  with  the  sixth  unknown 

7 


mass  fraction  determined  from  the  equation  /  Y,  a  l,  the  mass  fractions 


initially  differ  by  several  orders  of  magnitude  so  that  the  problem  was 
formulated  as  one  involving  six  differential  equations  and  the  summation 


condition  was  considered  a  check  in  the  accuracy  of  th©  results,* 


The  cases  studied  are  given  in  Table  II  along  with  the  mass  fraction 


Y4  in  each  case.  The  following  constants  w@r©  employed; 


*lt  was  found  that  die  summation  was  l  t  10"®,  Not#  that  if  the  element  mass 
fractions  are  considered,  only  four  differential  equations  need  be  considered. 
However,  die  same  objection  concerning  accuracy  applies. 


y  d  .v y  •  j  ;  A  .".-."A.  ,.v-'  :v  y  ■  v  .-r 


%  h 


V  ‘1-  - 


Aj  =18! 5.  23  cal/gm 
A^  =2687,9  cal/gm 
A3  =-2831.  3  cal/gm 
A^  =200.05  cal/gm 
A^  -  5938.  1  cal/gm 
A^  =  55,  873.0  cal/gm 
A_,  =911.78  cal/gm 
WI  =32 
W  =2 

W  =18 

3 

W  =28 

4 

W  5=I6 

w6=l 

W?  =  17 


=  .2  6403  cal/gm  °K 

Spz  -  3.  654  cal/gm  °K 

Cn  =.  56361  cal/gm  °K 
‘  3 

Cp  =.28382  cal/gm  °K 

CP5  cal/gm  °K 

c_  =4.  9  68  cal/gm  °K 
P6 

=  .43765  cal/gm  K 
Tr  =  1060°K 

Ru  =82.0  6  cc-atmos/mole°K 


& 

w 
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III  PRESENTATION  AND  DISCUSSION  OF  RESULTS 

The  results  of  the  cased  studied,  as  listed  in  Table  II,  are  presented  as  time 
histories  of  the  temperature  and  mass  fractions.  In  Figures  la  to  Id,  the 
temperature  histories  are  given  while  in  Figures  2a  to  2d  through  8a  to  8d,  the 
histories  of  the  mass  fractions  Yj  through  Yj  are  presented. 

The  history  of  the  reaction  can  be  characterized  by  three  periods.  Consider 
qualitatively  these  results;  first,  there  exists  an  induction  period  during  which 
there  are  significant  increases  in  the  mass  fractions  of  O,  H,  and  OH  not  observable 
from  the  graphs  but  given  by  the  computer  results.  At  the  end  of  this  induction 
period,  denoted  as  Tj,  there  occurs  a  second  period  during  which  a  rapid  change 
in  composition  occurs  with  relatively  little  change  in  static  temperature.  The  mass 
fractions  of  02,  H2  and  H20  come  close  to  their  equilibrium  values  during  this 
second  period.  However,  small  but  significant  amounts  of  atoms,  O  and  H,  and  of 
hydroxyl  radicals  OH  are  created;  these  amounts  are  in  excess  of  their  equilibrium 
values.  Since  the  enthalpy  per  unit  mass  associated  with  the  atomic  hydrogen  is 
from  20  to  200  times  larger  than  that  of  any  other  species,  this  is  by  far  the  most 
important  intermediate  in  the  system  under  consideration.  The  time  corresponding 
to  the  end  of  this  second  pei  iod  is  denoted  r2  «  During  the  final  and  most  extended 
period  of  the  reaction,  the  concentrations  of  the  intermediates  decay  to  their 
equilibrium  values  and  the  heat  release  associated  with  the  recombination  of  O  and 
H  and  the  decay  of  OH  leads  to  the  temperature  rise.  The  time  at  which  equilibrium 
is  essentially  achieved  will  be  denoted  r  .  This  general  scheme  for  the  reaction 
is  substantially  the  same  as  that  indicated  by  .Duff,  Schott  and  co-workers 
(References  2,  3  and  9)  for  the  reaction  profile  behind  a  detonation  wave. 
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Consider  more  quantitatively  the  numerical  results.  Given  in  Table  III  are 


the  times  i  =  1,  2,  3  for  the  various  cases;  these  were  estimated  by  examining 


the  time  history  of  the  mass  fraction  of  H.  The  time  t  was  selected  at  the  time 


when  an  observable  amount  of  H  appeared  on  the  graphs  (Figures  6a-d);  this  was 


Y^  «10  ,  The  time  T2  was  taken  to  correspond  to  the  peak  in  and  the  time  Ts> 

when  available,  was  determined  by  d  Y^/dt  — 0. 


The  ignition  delay  time  denoted  here  by  r  ^  and  predicted  by  the  approximate 


theoretical  analysis  of  Nicholls  (Reference  10)  has  been  interposed  between  the 


columns  for  Tx  and  T2  in  Table  III,  This  time  is  defined  as  that  corresponding  to 


an  inflection  point  in  the  curve  giving  the  time -history  of  Y  and  thus  should  be 

6 


between  Tj  and  r2.  The  equation  for  in  the  notation  and  with  the  rate  constants 


used  herein  is* 


rd=  Ru  V6Xl,oP<lol4>  e881°/T°  in(xi>0/X2 


+  jgn3  +  26340/T  (17) 


where  TQ  is  in  degrees  Kelvin,  p  is  in  atmospheres  and  Ru  =  82,06  cc-atm/mole-°K, 


Consideration  of  the  results  in  Table  III  indicates  that  within  o.  factor  of  two. 


Equation  (17)  predicts,  in  most  cases,  a  time  within  the  second  period  of  the  chemical 


reaction.  In  addition,  the  times  Tz  shows  a.  dependence  on  pressure  and  q  in 


satisfactory  agreement  with  either  Equation  (17)  or  Equation  (17a)**,  Thus,  the 


approximate  analysis  of  Nicholls  can  be  employed  to  give  the  dependence  on  pressure, 


temperature  and  oxygen  concentration  of  the  times  characterizing  the  early  phases 


*  Note  that  the  first  reaction  in  Equation  (1)  was  considered  by  Nicholls;  the  value 
given  by  Schott  (Reference  3)  for  the  rate  constant  was  used  in  Equation  (17).  It 
is  thus  implicitly  assumed  that  this  first  reaction  alters  only  slightly  the  time- 
history  of  the  composition, 

**  Note  that  Xj  q  =  (4.76  +  2 1?)-1 
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of  the  reaction  ;  indeed,  the  second  bracketed  term  in  Equation  (17)  is  roughly 

constant  for  the  initial  conditions  studied  here,  so  that  roughly 

„  8810/T 

T0e  ° 

T  ~ - -  (17a, 

d  Xl,OP 

There  seems  to  be  no  equivalent  analysis  for  the  estimation  of  the  overall 

time  t3  J  indeed,  one  of  the  desired  results  of  the  studies  of  chemical  models  would 

be  a  ready  estimate  of  this  overall  time,.  It  should  be  noted,  however,  that  Tj  is.. 

20  to  40  times  greater  than  tz  .  This  result  must  be  recognized  if  the  ignition  time 

delay  is  used  as  a  measure  of  the  chemical  time  as,  for  example,  in  Reference  6. 

The  extent  of  this  third  period,  during  which  the  heat  release  and  recombination 

occurs  implies  that  the  overall  reaction  time  is  sensitive  to  the  pressure,  roughly 
-2 

T3~  P  .  Thus  the  importance  of  pressure  on  the  overall  reaction  time  is 
emphasized,. 

There  are  shown  in  Figures  Ib-d  the  flame  temperatures  given  by  Reference  11: 
the  asymptotic  temperatures  given  by  numerical  integration  should  equal  these  flame 
temperatures.  It  will  be  seen  that  the  agreement  is  good,  so  that  the  assumptions 
with  respect  to  the  products  and  to  the  enthalpy-temperature  relations  may  be 
considered  satisfactory. 

It  is  perhaps  worth  noting  that  for  the  initial  conditions  corresponding  to  Case  7 
(p=  l  atmosphere,  Tq=  1600°K,  Tj~  1.8)  the  reaction  is  complete  in  roughly 
200  fj,  sec,;  thus,  if  a  mean  velocity  of  20,  000  ft/sec,  is  assumed,  a  combustion 
chamber  length  of  four  feet  would  be  required  to  achieve  complete  combustion. 
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IV  CONCLUDING  REMARKS 

A  study  has  been  made  of  the  reaction  history  of  hydrogen-air  mixtures  under 
adiabatic,  isobaric  conditions.  An  elaborate  chemistry  has  been  treated  herein  with 
a  view  toward  establishing  rather  accurate  reaction  histories,,  In  the  future,  simplified 
chemical  models  will  be  considered  and  the  reaction  histories  based  thereon  compared 
to  those  presented  here.  These  models  will  then  permit  approximate  analyses  of  the 
non-equilibrium  behavior  of  hydrogen-air  mixtures  under  flow  conditions  as,  for 
example,  in  exhaust  nozzles  and  mixing  regions. 

The  results  presented  for  a  variety  of  initial  temperatures,  pressures  and 
mixture  ratios  indicate  that,  during  the  early  phases  of  the  reaction  the  main  products 
achieve  close  to  their  equilibrium  values  with  little  change  in  temperature.  However, 
the  intermediates,  atomic  oxygen,  atomic  hydrogen,  and  the  hydroxyl  radical  exceed 
their  equilibrium  values.  The  remaining,  major  portion  of  the  reaction  time  involves 
the  decay  of  these  intermediates  and  the  chemical  energy  release  associated  therewith. 
Because  of  its  thermodynamic  importance,  i.e,,  its  large  enthalpy  per  unit  mass, 
atomic  hydrogen  appears  to  be  the  most  important  intermediate  for  the  reaction 
considered. 
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TABLE  I  -  RATE  CONSTANTS 


Reaction 


H  +  02 — 7*  OH  +  O 


3xlOwe 


14  -8810/ T 


O  +  H - »  OH  +•  H 

z 


3xl0He-4O3O/T 


OH  +  H2-^>H20  +  F 


3  x  10"  e 


14  -3020/T 


20H — r  H20  +  O 


Jx  1014  e~ 3020/T 


2H  +  M— >H  +  M 
z 


1  X  10 


H  +  OH  +  M— ^HzO  +  M 


1  x  10' 


H+O  +  M— *  OH  +  M. 


1  x  10 


20  +  M— >-0  +  M 


3x10 


OH  +  O  — »-H  +  O, 


2.48  x  10‘3e-66O/T 


OH  +  H  — ►  O  +  H, 


UxlO»e-i«WT 


H  O  +  H— >OH  +  H, 
2  i 


H20  +  O  — ^  20H 


1.33xlOl5e“l0>  950/T 
3.12*  l0l5e-l2,510/T 


H  +  M — ?  2H  +  M 
2 


H20  +  M — +  OH  +  M 


OH  +  M  — aH  +  O  +  M 


02  +  M  ->  20  +  M 


1.85  x  1019  „-54,  000/T 

e 

T 

9„66xl020  -62,  200/T 
- - -  0 

T 

8  x  1018  g-52,  500/T 

T 

2.  9  x  1019  -60,  600/T 


*  kg-kg  have  units  (moles/cc)"^  sec  * 


«1 

all  other  km  have  units  (moles/cc)  sec 


i\  >•  i  / 


<  "v\ 


V 
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TABLE  II 

1 

INITIAL  CONDITIONS  FOR  CASES  STUDIED 

■ 

Css® 

p,  atmospheres 

H 

o 

o 

pi 

V 

Y4  i 

i 

1 

0.  1 

1600 

0.2 

0.76249  j 

2 

0,  1 

1600 

Id  8 

0.72872  1 

3 

l 

1400 

0.2 

0.76249 

4 

1 

1400 

1.8 

0.72872 

5 

l 

1600 

0,2 

0. 76249 

6 

1 

1600 

1.0 

0. 74522 

7 

1 

1600 

1.8 

0.72872 

8 

l 

1115 

0.4 

0. 75810 

9 

1 

1155 

0.8 

0. 74947 

10 

1 

1187 

1.2 

0.74103 

11 

l 

1211 

1.6 

0.73278 

. 

! 

1  f 

I 
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TABLE  III 


TIMES  FOR  PERIOD  IN  REACTION  HISTORY 


Case 

T 

1 

T, 

d 

T2 

T 
'  3 

1 

10 

52 

27 

2 

10 

74 

52 

3 

3,4 

11 

7.0 

4 

5,0 

14 

12 

300 

5 

2,5 

5,2 

3.4 

6 

1.  6 

1.3 

3.8 

7 

1,8 

I.  1 

5,4 

200 

8 

25 

68 

39 

9 

22 

49 

32 

700 

10 

17 

44 

29 

11 

16 

41 

28 

AH  times  are  in.  microseconds 


APPENDIX  I 


EXPRESSIONS  FOR  7j  USED  IN  EQUATION  (4) 
h  s-kic6clfk9C7c5+k8c5c“k16clc 


K  ~“kZC5C2  +  k10°7c 6  "  k3c7c2  +  kllc3C6  +  k5C6c  “  k13C2C 


y  3  =  k  3C7C2  '  kllC3C6  +  V?  "  k12°3C5  +  k6C6C7C-*k14c3C 


y  =0 

4 


y5  =  k1c6ol  -  k9c7c5  -  k2c5c2  +  k10c7c6  +k4c7  -  k12c3c5 

2 

"  k7c6C5C  +  ki5c7c  "  2k8c5  c  +  2kl6clc 


^6  =  "klC6Cl  +  k9C7C5  +  k2C5C2  ”  k10c6c7  +  k3c7c2 


-klIC3C6  “  2k5C6  c  +  2kI3C2C  ’  k6c6c7c 


+  k14C3C  “k7c6c5c  +k15c7c 


r7  *V6C1  ~  k9°7C5  +k2C5C2  "  k10c7c6  ‘  k3V2 


+  k  c  c  -  2k  c  +2k,  c  c,  *  k  c 
U  36  47  12  35  6  6  * 


+  k14V+V6Vkl5V 
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FIG  f  T£  MPERATURE  HISTORY 


FIG.  I  TEMPERATURE  HISTORY 
c)  CASES  5-7 
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FIG..  3  HISTORY  OF  THE  MASS  FRACTION  OF  MOLECULAR  HYDROGEN 

b)CASE  3 


FIG.  3  HISTORY  OF  THE  MASS  FRACTION  OF  MOLECULAR  HYDROGEN 

b)  CASE  4 


FIG.  3  HISTORY  OF  THE  MASS  FRACTION  OF  MOLECULAR  HYDROGEN 

c)  CASES  5-7 


FIG.  3  HISTORY  OF  THE  MASS  FRACTION  OF  MOLECULAR  HYDROGEN 


FIG.  4  HISTORY  OF  THE  MASS  FRACTION  OF  WATER 

a)CASES  I  AND  Z 


FIG.  4  HISTORY  OF  THE  MASS  FRACTION  OF  WATER 
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c)CASES  5-7 


FIG.  4  HISTORY  OF  THE  MASS  FRACTION  OF  WATER 

d)  CASES  8-11 


FIG.  5  HISTORY  OF  THE  MASS  FRACTION  OF  ATOMIC  OXYGEN 

o)  CASES  I  AND  2 


FIG.  5  HISTORY  OF  THE  MASS  FRACTION  OF  ATOMIC  OXYGEN 

b)  CASES  3  AND  4 
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FIG.  6  HISTORY  OF  THE  MASS  FRACTION  OF  ATOMIC  HYDROGEN 

a) CASES  S  AND  2 


FIG.  6  HISTORY  OF  THE  MASS  FRACTION  OF  ATOMIC  HYDROGEN 

t?)CASES  3  AND  4 
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FIG.  7  HISTORY  OF  THE  MASS  FRACTION  OF  THE  HYDROXYL  RADICAL 
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